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Abstract 

The possible existence of a light CP-odd Higgs state in many new-physics models could lead to 
observable effects in the bottomonium sector. Experimental bounds on radiative T decays through 
such a pseudoscalar state and possible mixings with the states are reviewed. Combining these two 
effects, we set constraints on the properties of the CP-odd Higgs in the limit of small photon energy of 
BR(T — > 77V), that is on the pseudoscalar mass-range ~ 8 — 10 GeV. 

1 Introduction 

The bottomonium sector is well-known as a possibly sensitive probe of light CP-odd Higgs states ^4° [IH7] . 
Such pseudoscalar states have indeed vanishing V — V — A (V = Z, W) couplings, which allows them to 
circumvent most of the direct collider constraints (contrary to CP-even states). Several new-physics models 
provide a natural framework to embed these particles: MSSM with CP-violating Higgs sector 0[9], U(l)- 
extensions of the MSSM [TU], Little-Higgs models, (non-supersymmetric) Two-Higgs Doublet models [TTj . 
etc. 

As an illustrative case, we consider the example of the NMSSM (see [12J for a review), where a light CP- 
odd Higgs can emerge {e.g.) naturally from an approximate and spontaneously broken R or Peccei-Quinn 
symmetry. A phenomenological possibility in this model is the so-called "ideal Higgs scenario" , where the 
lightest CP-even Higgs h\ decays unconventionally into a pair of light CP-odd Higgs states with masses below 
the B — B threshold [13 20J so as to avoid LEP bounds on e + e~ — » Z + bb (since A Q — > BB is kinematically 
forbidden). This mechanism allows for comparatively light CP-even Higgs states > 86 GeV and can also 
lead to a successful interpretation of the 2.3 a excess observed in e + e~ — » Z + bb for M b i ~ 98 GeV [15] . The 
recent analysis of e + e~ — » Z + At in ALEPH data [2TJ constrains this scenario further but, according to [2"2"] . 
these bounds might also be circumvented by allowing for increased A — > cc decays; it is otherwise possible 
to simply increase the singlet component of the lightest CP-even state, which has reduced couplings to the 
SM-sector, e.g. to the Z boson. Note also that the presence of a light CP-odd state in the NMSSM spectrum 
is not necessarily associated with this "ideal Higgs scenario" and may as well accompany a heavier CP-even 
(hi) state (up to rrih 1 ~ 140 GeV in the NMSSM). Investigating the existence of such an A in low-energy 
signals is thus a necessary probe, complementary to Higgs searches at high energy. In this context, the 
bottomonium sector is of particular interest since light CP-odd Higgs states could be produced in T decays 
or mix with ^-states, provided their coupling to 6-quarks is sufficiently large [2"3"I|24) . 

Radiative T decays through a CP-odd Higgs state are constrained by several experimental limits, the 
most recent ones originating from T(1S) — > jA , A — » searches at Cleo III [55] and T(3S) — > 

7 A , A -> l+l~ searches at BaBar (I = [m, t). These experimental limits on a light NMSSM CP-odd 
Higgs have already been studied, in [51123] for the Cleo hi bounds and [22] for the BaBar bounds, and were 
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shown to constrain most of the significantly coupled region Xd > 0.5 — 1 (Xd being the reduced coupling of 
the A to down-type quarks). 

Such contraints should however be considered with caution in the region uia > 8 GeV, where BR(T — > 
7 A ) is not well controlled theoretically and could suffer from large corrections. Moreover, possibly relevant 
mixing effects with 775 (or xo) states could complicate the situation further 1,6,24 . Such limitations leave all 
the mass region wla ~ 8 — 10.5 GeV essentially unconstrained (or unreliably constrained) by these previous 
analyses and, in particular, the higher reach in mass of the BaBar bounds (m TT < 10.10 GeV) would seem 
impossible to exploit. 

Alternative searches for a light CP-odd Higgs signal through a breakdown of lepton universality in 
inclusive leptonic T decays have likewise failed so far [37]. Yet, for the same reasons as above, few conclusions 
can be drawn on the mass range rriA > 8 GeV. 

The mixing between the A a and ijb states was also studied and, in particular, it was suggested in |24j that 
this effect be responsible for a displacement of the observed rji,(lS) mass, which (depending on the QCD-based 
model) is generically lower than what was expected for the hyperfine splitting my (is) ~ m »;t,(is) [281125] . 

The purpose of this paper consists in investigating the impact of the mixing effect on radiative T decays 
with small photon energy, that is in the CP-odd mass-range wia ^8 — 10 GeV. We show that even when 
the direct decay T(15/35) — > jA°, which had been considered in [5,22,23 , is neglected, the additional and 
indirect contributions T(15/35) — > 7(77° ^4°) resulting from the mixing effect, lead to severe limits on 
Xd when compared to the experimental bounds from Cleo and BaBar. Given the large uncertainty on 
hadronic parameters and the limited validity of our approximations, such limits should however be considered 
cautiously and from a qualitative, rather than quantitative, viewpoint. In the following section, we summarize 
the status of BR(T —> jA ^ l + l~)) in view of the Cleo and BaBar limits, in the approximation of a 
pure A Higgs state. We will also present briefly the mixing effect between the A and the rjb states. In the 
third section, the consequences for BR(T — > 7t + t~) will be analysed in the limit where vtla ~ my, under 
the assumption of negligible direct decay T — > 7A . The numerical bounds on Xd will be discussed in the 
last section before a short conclusion. 

In the following, we will assume the existence of a NMSSM-like CP-odd Higgs state with mass itia < 10.5 
GeV and coupling to down-type quarks and leptons of the form "^^ d (where rrif is the fcrmion mass and 
v — (2\/2Gp)^ 1 ^ 2 is the electroweak vacuum expectation value), while the coupling to up- type quarks is given 
by ^J t& ^2 p ■ I n * ne NMSSM, Xd — cos 9 a tan/3, where cos 6 a quantifies the amount of doublet-component 
in A while tan/3 is the usual ratio of the doublet vacuum expectation values. Significant corrections to the 
Yukawa couplings are known to develop at large tan/3, due to loops of supersymmetric particles (see e.g. 
[3"U]). They would alter the relation Xd — cos 8 a tan f5 somewhat. Such effects will be neglected here, but 
we will see that constraints from the bottomonium sector can become relevant already at moderate tan/3. 

Moreover tan/3 > 1, so that the decay ^4° — » t + t~ is dominant (if kinematically allowed). In fact, the 
dominant branching ratios are essentially independent from tan j3, once tan j3 > 2 — 3 [22] . The results that we 
present in the following will assume the branching ratios at tan (3 — 5 (about 0.9 to 0.75 for BR(A° — > t + t _ ) 
in the range itia = 8—10 GeV; however, BR(A° — > t + t~) ~ 1 would have been a qualitatively acceptable 
approximation). Note that for tan/3 < 2, the coupling of the ^4° to 6-quarks is reduced and would generate 
little effect on the bottomonium sector (except for cos# 1 ~ 1). In the general case, the full (m \, tan 0) 
dependence can be retained. The branching ratio BR(A° —> r + r _ ) (also BR(A° —> /i + /i - )) that we use 
is obtained with the public code NMSSMTools [3T]. It was noted that the corresponding value can be 
slightly different from the output of Hdecay [22. . However, such minor effects will have little impact on 
the discussion that follows. Note also that the very-light mass region itia < 5 GeV is already essentially 
excluded by low-energy constraints [32 , such as B s — > or B — > X s /i + \i~ [33] . 

Beyond the NMSSM, this analysis should be essentially valid for any model with a light CP-odd state 
and based on a Two- Higgs Doublet model of type II (or even of type I), provided Xd is suitably chosen. 
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Figure 1: Constraints from BR(Y(nS) — > j(A° — > t + t~ / /i + on the mass / coupling plane of the A 
state. The curves correspond to upper bounds on Xd resulting from the experimental limits. The red, dark- 
and light-blue curves correspond, respectively, to Cleo limits on T(IS') — > i(A — > t + t~), BaBar limits 
on T(35) 7(A° r+r") and BaBar limits on T(35) ->■ j(A° m + aO- The correction factor F of eq. 
((T|) was assumed to be that shown in [23] and vanishes for ttia > 8.8 GeV. 

2 Direct bounds on BR{T -> 7 A°) and Mixing A - rj b 

In first approximation, Upsilon decays through a CP-odd Higgs are described by the so-called Wilczek 
formula |34) : 

BR(T(nS) -» 7 A°) _ Gf-mgXj / _n 
BR(Y(nS) -> ~ v^ra ^ ~ m| 

This formula was first established in the non-relativistic approximation and the hard-photon limit, i.e. 
rriA <C mi' Several corrections, due to relativistic, QCD and bound-state effects, were then included within 
the factor F (see [35] for a summary). Yet our control over this factor fails for tua > 8 GeV, in the limit of 
low-energy photons, where large corrections from e.g. bound-state or relativistic effects are expected. Within 
the known approximate computations, the correction factor F vanishes for ttia > 8.8 GeV (depending on 
the 6-mass), hence leaving this mass range apparently unconstrained |23j . One can indeed expect the A 
contribution to radiative T decays to vanish in this limit, since, for large wavelengths, the photon simply 
probes an overall neutral state. Yet, the cancellation for rriA > 8.8 GeV overestimates probably this feature. 
In fact, we will show in the following sections that the mixing effect, for itia ~ w^o, leads to sizable 
contributions, which are not captured by eq. (|TJ) . 

Comparing eq. ([T]) with the experimental limits from Cleo [25] and BaBar (53] on BR(T(1S/3S) — > 
j(A° — > ii + \i~ /t + t~)), one obtains upper bounds on Xd as a function of ttia [5l l22l[23j . as shown in fig. [TJ 
For the mass region beneath ttia ~ 8.8 GeV, only moderate to small values of Xd < 0.5 are allowed. Note 
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Figure 2: Reduced coupling of the light CP-odd Higgs, Xd = cos 9 a tan/3, as a function of tan/3, in units of 
-p^-, for the R (red curve) and PQ (green) limits. A and (i Bft are NMSSM parameters [12] . with v ~ 175 GeV, 
A < 0.8 and > 100 GeV. We thus obtain that Xd is generically small in these limits. 



that this region of weak coupling is however that which is usually favoured from the theoretical point of 
view: the R and Peccei-Quinn symmetry limits of the NMSSM would naturally predict small Xd (see fig. 
[2]). On the other hand, such weakly-coupled A will prove difficult to observe experimentally. 

As we discussed in the introduction, the mass range itla ~ 9 — 10.5 GeV is (due to the smallness of the 
factor F in eq. ([IJ, within the known computations) free of the latter constraints, so that large values of Xd 
seem a priori allowed. In this range however, one expects to find the rjb(nS), n = 1,2, 3, states. Such states 
carry the same quantum numbers as the CP-odd Higgs, so that a mixing A — ijb develops. This scenario 
was studied in [l][6l[24]. We shall use the notations of [24]. The mixing can be described by an effective 
mass-matrix for the rjb — A states: 
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(3) 



The quoted value for m^o^g) is the perturbative QCD (pQCD) prediction. Slightly different predictions 
can be obtained in other approaches [29]. The uncertainty on the heavier 77° masses will be neglected. The 
off-diagonal terms can be computed in a non-relativistic quark-potential model [T]|ni[23] : 
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Figure 3: Relations and constraints in the plane (mj,!^ originating from the splitting effect in the pseu- 
doscalar mass matrix and the observed 777, mass. The red curves represent upper bounds on Xd'. the full line 
corresponds to the conservative limits —30 MeV < m^o( lg j — m oba < 54 MeV; the dotted one corresponds to 
the pQCD 2cr range. In the domain between the orange dashed curves, the splitting between the observed 
rib mass and the pQCD prediction is generated by the mixing effect within f a. The green curve is the 
"favoured" line obtained with the central value of the pQCD prediction 24 . 



where the r^inS) wave-functions R Vb ( n s) were estimated through those of the Y(nS) states, which can be 
extracted, in turn, from the leptonic decays T(nS) —> [36]. The other off-diagonal entries vanish as a 
result of the orthogonality between 77° states. The physical mass-states are denoted as: 



m = Pn r%(lS) + Pr2 %°(2S) + P l3 77° (3S) + P l4 A 



(5) 



The only experimental information at our disposal concerning the mass matrix of eq. ^ originates from 
the observation at BaBar, in T(3S) [37] and T(25) [35] decays, of a state with mass m obs = 9390.9±3.1 MeV 
[38] which was interpreted as the 77° (IS). It is remarkable that this value for the mass is in slight tension with 
the pQCD prediction (eq. ([3])). Within most of the QCD models, one expects a slightly smaller hyperfme 
splitting [5S] (although possibly within la). As a consequence, the mixing effect with the A was suggested 
as a possible interpretation of a displaced 77° (IS) mass [53]. The requirement that mf bH be an eigenvalue of 
the mass matrix @ leads to the relation: 



5 TOj 



X(3S) 



(6) 



For given QCD-predicted r\b masses (see eq. ©) and the mixing elements of eq. (@|, eq. ([5]) determines rriA 
in terms of Xd- Varying the input from eq. (|3J) and (jU within the error bars, this relation results in bounds on 
the plane (jua, Xd) [23) . These limits originate simply from the observation that, in the presence of a mixing, 
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Figure 4: Contributions to the radiative Upsilon decay to a Higgs pseudoscalar. 

a) Direct decay, corresponding to the description of eq. (JXJ) - 

b) Additional mixing contribution due to the A — rjb mixing. 

a mass-splitting is generated so that the observed mass cannot coincide with a diagonal entry of the mass 
matrix (to which it couples). Assuming the pQCD result to be valid, one can define a favoured region in the 
plane (m^, Xj) where the observed mass can be reproduced within 1 a. This is illustrated in fig. [3] The 2 a 
limits are also shown. However, given the larger number of predictions in alternative QCD-models, we retain, 
in what follows, the bounds resulting from the more conservative range —30 MeV < rn^o^g^ —m obs < 54 MeV. 



3 Radiative T decays in the mixing scenario 

B R(T (nS) — > 7A ), as we presented it in the previous section, corresponds to the diagram of fig. 3^,. In 
the presence of a ^4° — rjt, mixing, however, a second contribution, as in fig. \^p, arises. This contribution 
was already described in [B] under the denomination "resonant decay" . Assuming that the contribution of 
fig. @Ji vanishes for small A°-T mass differences, as we discussed in section^ the diagram of fig. HJd is then 
dominant and could lead to a significant effect. In this section, we will explicitely neglect this diagram of 
fig. [4^, when ttla > 8 GeV and focus exclusively on the contribution b due to mixing. 

Let us first consider the transition between pure bottomonium states. In the non-relativistic approach 
to quarkonium bound states, BR(T(nS) — » Jil®(jS)) can be written as (3"9"] : 

BR(T(nS) -+ 7% °(iS)) = - \ a / b \I nj \ 2 k 3 ■ 0(fc) (7) 

where a — (137. 036)^ is the fine-structure constant, eb = —1/3, the &-quark charge, k, the photon energy 
and r T („5), the Y(nS) total width (9 is the Heaviside distribution). We shall use r X (is) = (54.02±1.25) keV 
and Tx(3S) = (20.32 ± 1.85) keV [35] ■ The transition factor I n j is defined as: 

I nj = (v° b US)\jo(kr/2)\?(nS)) (8) 

where ja(x) = sa*. Taking tq ~ 1.4 GeV" 1 as a typical scale of confinement, we see that this function can 
be safely expanded in the limit of low-energy photons as jo(fcr/2) ~ 1 — (fcr/2) 2 /6 + . . .. Using that the wave 
functions for the rfl and T states are almost identical, one obtains that in first approximation I n j ~ 6 n j: 
nS — > nS transitions are a priori favoured. The resulting k 3 dependence in eq. (0 is then characteristic of 
point-like interactions of the bottomonium states, which could be anticipated for large wavelengths. 

Note however that for QCD bb states, the photon energies are much smaller for a nS — > nS transition 
than for a 35 — > 1S/2S transition, resulting in comparable branching ratios. The coefficients I n j, n =/= j, 
were estimated in QCD (see e.g. [3S]), taking into account relativistic corrections, but such predictions range 
over one to two orders of magnitude. Fortunately, this situation will be slightly changed in the case of mixed 
states (see eq. (fTQ)) below) and it will be possible to focus on the dominant nS — > nS channel. For the time 
being, let us use that I n j — S n j and that photon energies are small, and parametrize: 

I n j — S n j k • R n j (9) 
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So far the R\^ coefficients are simply unknowns, which might be extracted from the QCD computations of 
the I n j. Under the assumption that the nS — > nS transition dominates, it will be possible to neglect the 
R 2 nj , n^j. 

In fact, for large mixing (that is, for tua in the vicinity of the 77° masses), processes involving the CP-odd 
Higgs or the rfl states are more adequately described in the context of the mixing formalism which we have 
presented in section [5] In particular, the production of any mass state rji, i = 1 — 4, in radiative T decays 
may be obtained (under the assumption that the contribution of fig. 2^, is negligible) by a generalization of 
eq. © as: 

A 2 3 
BR(T(nS) -+ lm ) = - " 6fc kf ■ e(h) Y, Zjlm PijPil (10) 

2 _ 2 

Note that the photon energy is hi = — ''' f° r ah the terms on the right-hand side of eq. (TlOl) . 
including pure nS — > nS transitions as well as nS — > nS, n ^ h, transitions or interference terms. Under 
such conditions, the nS — > nS transition dominates eq. (|10[) due to its larger transition factor. In the 
following, we will neglect the interference terms in eq. (fTUJ): such terms cannot be large unless the rji has 
simultaneously large r)®(jS) and ri®(lS), j 7^ I, components, which (almost) never happens since the 77° states 
are orthogonal; moreover, both transition factors I n j, I n i, j 7^ I, would have to be large, which contradicts 
eq. © for small photon energy. Eq. ((T0|) thus reduces to: 

3 

BR(T(nS) -> im ) ~ J2 P l BRi(T(nS) -> j V ° b (jS)) (11) 
i=i 

where BRi(T(nS) —> "fT)b(jS)) is given by eq. ([7]), while the index i reminds us however that the relevant 
mass for the rf state is that of the mass state r\i . Note that the simplification from eq. (|TU)) to eq. (jlip will 
work well due to the dominant nS — > nS channel. This also implies, however, that we restrict ourselves to 
small photon energies (hi), so that I n j ~ 8 n j remains a good approximation: therefore, we will require in 
the following that ki < 1 GeV. 

Then, the decay of the mass state rji to a tauonic pair can be assumed to proceed through the A 
component |24j : 

BR(th -> t+t-) = P ^( A °^r+r-) 

where the widths r,^^ 2,3s) can be estimated to 10, 5 and 5 MeV respectiveljQ- T A o and Ti(A° r + r _ ) 
can be computed. Again, the index i indicates that the partial width should be calculated using the mass of 
the state rji instead of that of the diagonal A . The BR(rji —> r + r _ ) has already been proposed as a probe 
for light pseudoscalars in [4"T] . 

We can eventually write the radiative T decay to a tauonic pair as: 

4 

BR(T(nS) ->7tV) ~ J2 BR ( T ( nS ) ~> Th) ' BP (Vt t + t-) (13) 

»=i 

where BR(T(nS) — > 777^) and BR(rji — > r+r - ) are given respectively by eq. (TTT|) and (TT^j) . We stress that, 
in this analysis, the main contribution to eq. ([13")) is associated with the subprocess: 

T(nS) -> 7 ( 7]b(nS) — A -> r+r" ) (14) 

the nS — ^ n5 transition being favoured by its large transition factor while the mixing effect allows it to 
access photon energies comparable to those of the other transitions. 

At this point, neglecting the direct contribution T — > 7A of fig. 0^, was the most important assumption. 
If the factor F of eq. ((T|) is indeed small in the considered mass-range (tua ~ 8 — 10 GeV), as its cancellation 
in known computations seems to indicate, we may infer that the amplitude of fig. 2£i remains small, so that 



1 After r^o (nS) /r^ (nS) ~ ( mb /m c )[ as (m b )/a s (m c )] 5 ~ 0.25 - 0.75 [40]. 
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our assumption in neglecting the direct contribution is justified. In practice, the limits on Xd that we obtain 
(in section using the indirect decay (fig. |Dd) only, are comparable in the mass range ttia = 8—10 GeV 
to those that eq. ([T]) alone would have given with F ~ 0.5 — 1 (for the BaBar bounds) or F ~ 1 (for the 
Cleo bounds): therefore, we conclude that a factor F < 0.1 should be a sufficient condition to ensure that 
the contribution of fig. 0}d captures the dominant effect, hence that our approach is valid. Note that, among 
the neglected terms, the non-interferent contribution given by eq. (JlJ could be directly included in eq. (jlip 
(if F were satisfactorily kwown) and, being positive, would only strengthen the bounds that we obtain in 
section @] However, if the amplitude of fig. 0Ji were non-negligible, the interference between the diagrams 
of fig. 2k and |3Jd could be significant: since its sign is not controlled, this would spoil the validity of our 
analysis. 

A remark concerning the choice of the diagonal r]®(lS) mass is necessary. In the presence of the mixing 
structure of eq. ([2]) and knowing one of the eigenvalues m obB , one cannot choose rn v o^ 1S y ttla and Xd 
independently (or m obs will not be recovered in the general case). Therefore, we will determine m^o^g) in 
terms of rriA, Xd and m obB through eq. (JB]). Should m^o^g) reach unrealistic values (that is, if the condition 
—30 MeV < m^o^s) — m oba < 54 MeV is violated), the corresponding points will be excluded through the 
mixing contraint (see end of section [5] and fig. [3]). 

Finally, the advantage of considering the three rf^ states simultaneously lies in a more accurate and 
conservative description of the T(nS) — » 777^ transition: whereas the 77^ (j — nS) component is the most 
relevant one, the presence of the other 77° (j ^ nS) components reduces the size of the corresponding mixing 
element, hence alleviates slightly the impact of the experimental bounds. 

We thus conclude that, if only through the contribution of fig. HJd, the experimental bounds on BR(T(nS) —> 
jt + t~) should translate into constraints on the plane (mA,Xd), even in the limit uia ~ 9 — 10 GeV. 

4 Bounds on the light CP-odd Higgs 

In this section, we present the consequences of the analysis of section |3J for the light CP-odd Higgs in the 
limit of small photon energies. Given the large uncertainties on the transition factors I n j, one cannot expect 
to obtain strict bounds. Instead, we evaluate eq. (|13[) under several approximations and compare the results 
to the experimental limits. 

The BaBar limits on BR(T(3S) -> 7r+r~) have the highest reach in mass (m TT < 10.10 GeV, limited 
kinematically by the T(35 l ) mass and the requirement of observable, not-too-soft photons), hence are of 
particular interest. As a first approach, we will consider only the 35 — > 35 transition with the approximation 
that I ^33 1 2 — 1 , neglecting the coefficients i? 2 ^ . Then we will continue to neglect the transition factors I n j , n 7^ 
j, but assume | -Z33 1 2 is reduced down to 1/2. Finally, we will use eq. © with the values - 0.02 GcV~ 2 , 
R\ 2 ~ 0.25 GeV -2 and i? 2 3 ~ 0.02 GeV -2 . This choice of R\ x allows to reproduce approximately the 
measured BR(T(3S) -> jrj obs ) = (4.8 ± 1.3) ■ 10~ 4 [37] in the limit where no mixing is assumed. i?§ 2 
was chosen so as to have the same approximate ratio with i? 2 ^ as in |39] . i?§ 3 is simply the order 1 term 
in the expansion of j'o(Axo/2), ro being the typical scale of confinement. The choice of these values is 
of course arguable: the main purpose here consists in checking that the leading effect is captured by the 
35 — > 35 transition. The corresponding limits on the (m,A,Xd) plane are shown in fig. [5ji. We varied all 
the experimental input within 2 a before extracting the bounds. The 77^ widths were however kept at 10, 5 
and 5 MeV respectively: their effect is studied separately, in fig. [5}d. We observe that values of Xd > 1 — 2 
are constrained severely and more or less equivalently in the three cases we considered and we regard this 
fact as a sign of robustness of the corresponding limits: despite significant variations of the I n j factors in the 
three cases under consideration, the excluded region in the (m^.Xj) plane remains essentially unchanged. 
The underlying reason for this is related to the dominant 35 —> 35 transition. Fig. (5}d leads to a similar 
observation: increased 77° widths alleviate slightly the bounds as the total widths of the 77, state is increased 
(leading to a smaller branching ratio into t + t~). Yet despite a doubling of the (conjectured) 77" widths, the 
bounds on the properties of the A remain little affected. The dip for tua ~ ™ obs corresponds to unacceptable 
values of m^o^gy For tua < 9.3 GeV, the lightest state (dominantly A ) leads to larger photon energies so 
that the plotted bounds become unreliable. 

Although our approach to the BaBar limits should fail for 771^4 < 9.3 GeV, the Cleo limits can then be 
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Figure 5: Upper bounds on Xd as a function of uia due to the mixing effect in BR(Y(3S) — > jt + t~) (using 
BaBar data). 

a) Influence of the I n j coefficients on the bounds. ^ri {i,2,3S) — 10' 5, 5 MeV respectively. The full light-blue 
curve assumes i? 2 „ = and I/33I 2 = 1/2; the dashed middle-blue curve corresponds to i? 2 „- = and I/33I 2 = 1; 



the dotted dark-blue curve assumes the R\^ as described in the main text 

and |/ 33 | 2 = 1, wiili V 



b) Influence of the T)® widths on the bounds. J2| - 
(middle-blue curve) and f^o^ ,2,3s) = 20, 10, 10 MeV (pink curve). 



■«7?(1,2,3S) 
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Figure 6: Upper bounds on Xd as a function of tha due to the mixing effect in BR(Y(IS) — > 7t + t~) (using 
Cleo data). 

a) Influence of the I n j coefficients on the bounds, r^o^ ,2,3s) — 10, 5, 5 MeV respectively. The full orange 
curve assumes i? 2 - = and |/n| 2 = 1/2; the dashed red curve corresponds to R\^ = and |/n| 2 = 1; the 
dotted dark-red curve assumes the R\^ as described in the main text. 

b) Influence of the 77" widths on the bounds. R\^ = and |/n| 2 = 1, with ^0^,2,35) = 10,5,5 MeV (red 
curve) and r^o( 12 ,3g) = 20, 10, 10 MeV (brown curve). 
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Figure 7: Conclusion: upper bounds on Xd for a pseudoscalar A in the mass-range 8 — 10.1 GeV. 

- in red, from BaBar, resulting from the direct T decay to a ^4° state (fig. [3^,): similar to fig. [TJ 

- in green, from Cleo, resulting from the mixed contribution (fig. HJd); the full line corresponds to the 
approximate validity range, while the dotted curve extends to regions where the condition on the photon 
energy k < 1 GeV is not satisfied (i.e., where the assumption I nn ~ 1, hence the corresponding bounds, are 
unreliable) ; 

- in blue, from BaBar, resulting from the mixed contribution (fig. Eb); the full line corresponds to the 
approximate validity range, while the dotted curve extends to regions where the condition on the photon 
energy k < 1 GeV is not satisfied. 



applied. The favoured transition is now 15 — > 15. As before, we consider the cases \hi\ 2 — 1, \In\ 2 — 1/2 
and finally R\ x ~ i?§ 3 , R\ 2 ~ ^321 ^13 ~ ^31 • The result is shown in fig. [5] and constrains Xd < 2 — 4 in 
the range tua ~ 8.4 — 9.4 GeV. Again, we observe little variations of the general aspect of the bounds in the 
several considered cases. 

We conclude this analysis with fig. [7]where we collect all the previously discussed limits on the (tua, Xd) 
plane in their approximate range of validity. (We chose the case I nn = 1, I n ^j = 0.) We stress that, 
due to the limited control on the coefficients I n j (and on the A — rjb interference), the precise limits on 
this (rriA,Xd) plane should be considered with caution, that is, as a qualitative trend, rather than a strict 
exclusion boundary. In view of the results of fig. [5] and [6] however, we conclude that values of Xd beyond 
~ 2—3 should lead to unacceptably large BR(T(nS) — > 7 t + t~), at least within the range of the experimental 
constraints. 

Concerning the possibility of explaining the tension between the pQCD prediction and the observed 
%(15) mass, most of the "favoured" region, where this discrepancy can be interpreted (at least within la) 
by the presence of the A , is now excluded. The variations of BR(T(3S) — > 7t + t~) in the approximation of 
eq. (|13[) . along this "favoured region", are shown in fig. [5] (with Xd now determined in terms of tua through 
cq. ([5])). Only two mass ranges remain, where the tension between the observed and the pQCD-predicted 77° 
masses may be interpreted as an A effect: in the first case, tua ~ 9.39 GeV, the required Xd is sufficiently 
small to escape the BaBar bounds. This coincides with the region of natural NMSSM couplings. Note 
however that the acceptable mass range is now only of a few 10 MeV around m obs , which supposes a very 
fine accidental relation between the A and 77° (15) masses. The second region tua > 10.1 GeV is beyond 
the reach of the BaBar searches and requires large Xd ~ 20. We must therefore regard the possibility of 
generating the 1% mass shift in the context of a light CP-odd Higgs only as a marginal possibility. 

To summarize, we reviewed experimental constraints from Cleo and BaBar on radiative T decays 
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Figure 8: Mixed — A contribution to BR(T(3S) — > jt + t~) compared to the BaBar limits (dark curve) 
in the "favoured" scenario of (2^- The result is shown along the "favoured" line (green curve), and the two 
1 a boundaries (orange curves). Except for a few 10 MeV around m oba ~ 9.39 GeV and the large mass region 
(tua > 10.1 GeV), unconstrained by BaBar, the essential part of the mass range in this scenario would 
lead to excessive contributions. 

through a pseudoscalar. Even though the pure Higgs contribution is little controlled in the limit of low- 
energy photons, we showed that the mixing effect should lead to a sizable contribution, which can be used 
as a theoretical estimate for BR(T(3S) — > 7t + t"). Taking this estimate seriously, we obtain that the 
reduced coupling of the ^4° to b quarks must be small, Xd < 2, even in the range tua ~ 9 — 10.1 GeV. Note 
however that this region with small couplings is that which is the most naturally achieved in the NMSSM (in 
e.g. the R or Peccei-Quinn symmetry limits). The possibility of interpreting the small tension between the 
observed and pQCD-predicted 77° (15) masses through a mixing with the A is nevertheless largely narrowed 
and confined to the ranges wla ~ m obs or > 10.1 GeV. 
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